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A  BUMPY  TORUS  BETATRON 


I.  INTRODUCTION 

Conventional  betatrons1  are  current-limited  due  to  the  defocusing  effects  of  space  charge  at  injec¬ 
tion.  In  recent  years  there  have  been  several  renewed  attempts  at  overcoming  this  (rather  severe) 
space  charge  limit.  Specifically,  there  have  been  high  current  conventional  betatrons  proposed2  which 
employ  high-energy  injectors  as  well  as  so  called  modified  betatrons3,4  which  employ  a  toroidal  mag¬ 
netic  field  to  prevent  space  charge  blow-up  of  the  beam.  In  both  of  these  cases  however,  a  mismatch 
between  the  injection  energy  and  vertical  field  of  a  few  percent  or  so  will  cause  the  beam  to  hit  the 
wall,  a  matter  of  some  concern  in  a  high  current  device.  The  maximum  allowed  error  in  the  vertical 
field  is  typically  on  the  order  of  a  few  gauss  in  designs  which  have  been  considered.  Recently  it  was 
shown3  that  the  combination  of  an  /  -  2  stellarator  field  and  ordinary  weak  focusing  betatron  field 
results  in  a  strong  focusing  high-current  betatron  or,  *stellatron,"  with  a  large  energy  bandwidth.  Such  a 
configuration  has  the  advantages  of  relaxing  the  vertical  field  and  injector  tolerances.  In  addition,  the 
strong  focusing  introduces  a  threshold  for  the  negative  mass  instability,  so  that  this  instability  does  not 
operate  at  injection  (though  other  fast  growing  resistive  or  kink  modes  may  occur  below  the  negative 
mass  threshold).  In  this  note  we  report  analytical  and  numerical  results  on  the  bandwidth  and  stability 
of  an  alternative  strong-focusing  scheme,  namely,  a  combination  "bumpy  torus"  and  betatron  field, 
corresponding  to  the  /  —  0  stellatron. 

II.  DISCUSSION 


The  bumpy-torus  betatron  field  consists  of  a  superposition  of  an  /  -  0  stellarator  field  and  the 
field  of  a  conventional  betatron.  Near  the  minor  axis  at  r  -  r0,  z  -  0,  this  field  has  the  form 

B,  —  —  nyBm  +  y  8 B,  mx  sin  mO 

ID 

Bt  “  Bt0  (1  +  cos  m9)  (1) 

Bio 

B,  —  B„  (1  -  nx)  +  y  8 B$  my  sin  mB 

where  x  =  (r  —  r0)/r0,  y  =  z/r„,  9  is  the  azimuthal  angle,  n  is  the  vertical  field  index,  and  m  is  the 
number  of  bumpy-torus  field  periods  around  the  torus.  B„,  Bta,  and  bB$  are  constants. 


Treating  the  self  fields  of  the  beam  by  a  simple  cylindrical  model,  we  find  the  equation  of  motion 
for  a  test  particle  within  the  beam  is,  in  the  paraxial  approximation,  for  n  —  1/2, 


~~-f  +  — 7  12  —  4a,  +  82(1  +  e  cos  2  ffm)2]  $ 
d9m  m 
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where  9m  as  m9/2,  tjr  s  Or  +  (y)  exp  l(ib/2m)  (2 9m  +  <  sin  2  •„,)],  b  =  BtJB„ ,  c  =  8 BJBt0,  P„ 
is  the  momentum  of  a  particle  which  would  circulate  on  the  minor  axis,  bP  is  the  "momentum  error," 
n,  =  <u|/(2y£n£)  where  ub,  Cl„  are  the  beam  plasma  frequency  and  the  vertical  field  cyclotron  fre¬ 
quency,  respectively,  and  y0  —  (1  +  {Pjmc)1)111.  We  are  interested  both  in  the  solution  to  the  homo¬ 
geneous  part  of  Eq.  (2),  which  will  give  orbital  stability  criteria,  as  well  as  in  the  solution  to  the  inho¬ 
mogeneous  problem,  which  will  give  the  momentum  compaction  of  the  machine. 


The  quantity  n,  appearing  in  Eq.  (2)  describes  the  (net  defocusing)  effect  of  the  self  electric  and 
magnetic  forces  of  the  beam.  Since  it  depends  on  beam  density  and  therefore  on  the  beam  minor 
radius,  n,  will  in  general  vary  with  azimuthal  angle  9  around  the  device  in  a  manner  governed  by  the 
standard  beam  envelop  equation.  Consequently,  when  the  beam  envelop  is  stable,  we  expect  n,  to 
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behave  as  rt,(9)  ~  n„  +  nsl  cos  mff  +  ...  but  we  shall  assume  here,  for  simplicity,  that  2 eb2  »  nsi 
so  that,  in  Eq.  (2),  n,  may  be  adequately  approximated  by  its  average  value. 


f 


k 


Equation  (2)  is  a  Hill  equation,  which  has  characteristic  bands  of  stability,  as  shown  in  Fig.  1. 
The  boundaries  of  the  stable  regions  have  been  obtained  numerically,  using  standard  methods.6  The 
shaded  regions  in  the  figure  are  unstable  portions  of  the  plane,  e  vs  b/m,  for  the  case  n,  -  30  and 
m  -  30.  The  intersections  of  the  unstable  regions  with  the  abscissa  are  given  by 

(b2  +  2—  4na)/m2  -  q2  where  q  -  0,  1,  2  ... 

which  is  the  condition  that  the  transverse  rotation  frequency  of  particle  within  the  beam  is  an  integer 
multiple  of  the  "focusing  frequency,"  m  Cl  „  —  a  condition  which  allows  resonant  transfer  of  energy 
from  the  longitudinal  to  transverse  degrees  of  freedom  and,  consequently,  exponential  growth  of  the 
betatron  oscillation  amplitude. 


As  B2  is  increased  during  acceleration,  one  typically  would  not  wish  to  increase  Bt  simultaneously 
since  this  would  require  significant  additional  energy  storage.  The  result  is  that  the  operating  point  of 
the  accelerator  will  move  from  right  to  left  in  Fig.  1.  Consequently,  the  accelerator  should  be  run  in 
the  left-most  stable  band  to  avoid  crossing  unstable  bands.  These  considerations  require  m  >  b  at 
injection  and  force  the  use  of  a  large  number  of  field  periods  in  the  design  of  the  strong-focusing  sys¬ 
tem.  The  left-most  unstable  band,  corresponding  to  q  —  0,  is  due  to  the  beam  space-charge  and  rapidly 
disappears  during  acceleration  since  the  self-field  index,  n„  is  proportional  to  y? J,  where  y0  is  the  rela¬ 
tivistic  factor.  The  left-most  stable  band,  therefore,  becomes  broader  during  acceleration;  the  first 
stable  band  is  at  its  most  narrow  at  injection,  when  y0  is  smallest. 


We  next  consider  the  important  question  of  containment  of  particles  whose  average  momentum  is 
not  matched  to  the  vertical  betatron  field,  i.e.,  the  question  of  the  momentum  compaction  of  this 
configuration.  In  order  to  address  this  question  we  have  examined  numerically  the  behavior  of  single 
particle  orbits,  neglecting  beam  self  fields  but  employing  the  full  Bessel  function  representation  of  the 
I  -  0  focusing  field.  Figure  2  shows  the  allowed  mismatch,  BP/P„,  plotted  against  e  =  8 BJB,  for 

Bt0  -■  2 kG,  B„  -  118G,  n  -  y,  r0  -  100  cm  and  m  -  30.  This  plot  is  generated  numerically  by 

launching  particles  on  the  minor  axis  along  the  toroidal  direction  with  various  amounts  of  mismatch. 
The  figure  shows  the  largest  mismatch  for  which  the  calculated  orbits  are  contained  in  a  10  cm  minor 
radius  chamber.  Containment  of  particles  with  a  mismatch  of  ±  20%  is  obtained  for  c  =  0.2.  We 
stress  that  the  momentum  compaction  of  this  configuration  is  due  to  the  alternating  gradient  field  of  the 
"bumps,”  though  the  phase  shift  per  period  is  dominated  by  the  average  value  of  the  toroidal  field. 
Using  Eq.  (2),  with  n,  —  o,  a  perturbative  calculation  valid  for  small  values  of  c,  of  the  momentum 
compaction  factor,  gives 


8  r/rt 
8  P/P0 


(3) 


which  holds  only  for  m  >  b.  One  sees  in  Eq.  (3)  the  helpful  effect  of  a  bumpy  torus  field. 


In  conventional  betatrons,  resonances  are  automatically  avoided  by  increasing  the  particle  momen¬ 
tum  and  the  vertical  magnetic  field  in  synchronism.  The  introduction  of  non-synchronous  fields  (a 
fixed  toroidal  field,  for  example)  makes  the  betatron  wavelengths  energy  dependent,  which  can  lead  to 
the  crossing  of  resonances  driven  by  field  errors  during  acceleration.  As  in  all  strong-focusing  devices, 
the  occurrence  of  orbital  resonances  plays  an  important  role  in  the  operation  of  the  bumpy-torus  beta¬ 
tron.  Using  the  Floquet  solutions  to  Eq.  (2)  it  is  possible  to  obtain  a  condition  for  the  integer  reso¬ 
nances,  when  space-charge  effects  may  be  neglected: 


iir)  •  cos 


8  +  2*1 
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m 
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where  is  the  solution  to  Eq.  (2)  with  BP  a  0  satisfying  ifn(0)  a  1,  ^rj(O)  -  0  and  where  k  is  an 


integer,  the  Fourier  component  number  of  the  dipole  field  error.  Equation  (4)  provides  the  basis  for 
numerical  calculation  of  contours  in  the  stability  plane  on  which  Eq.  (4)  is  satisfied  for  a  given  Ar,  an 
example  is  given  in  Fig.  3. 

If  all  the  fields  cannot  be  made  synchronous  with  the  particle  energy,  the  effect  of  resonant  insta¬ 
bilities  might  be  minimized  by  making  the  energy  gain  per  pass  large.  Other  possibilities  for  coping 
with  resonance  crossings  are  currently  under  investigation. 

III.  SUMMARY 

In  conclusion,  we  find  the  spatially  alternating  transverse  magnetic  field  gradient  associated  with  a 
bumpy-torus  leads  to  a  potentially  interesting  strongly-focused  accelerator  configuration  which  is  seen  to 
have  a  region  of  stable  orbits,  and  to  have  a  significant  bandwidth  in  allowed  mismatch  between  the 
vertical  magnetic  field  and  the  particle  momentum. 
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